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Microwave-induced excess quasiparticles in superconducting resonators 
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We have measured the number of quasiparticles and their hfetime in aluminium superconducting microwave 
resonators. The number of excess quasiparticles below 160 mK decreases from 72 to 17 //m~^ with a 6 
dB decrease of the microwave power. The quasiparticle lifetime increases accordingly from 1.4 to 3.5 ms. 
These properties of the superconductor were measured through the spectrum of correlated fluctuations in 
the quasiparticle system and condensate of the superconductor, which show up in the resonator amplitude 
and phase respectively. Because uncorrelated noise sources vanish, fluctuations in the superconductor can be 
studied with a sensitivity close to the vacuum noise. 



The promise of a long quasiparticle lifetime and a long 
coherence time makes superconducting circuits popular 
for use in radiation detection and quantum computa- 
tion. At low temperature the number of quasiparticles 
in a superconductor should decrease exponentially. Ex- 
cess quasiparticles were recently suggested to limit the 
coherence time of superconducting qubits^ and the tun- 
nelling rate in single-electron transistors^,. Recently, the 
quasiparticle lifetime in a high-quality aluminium super- 
conducting resonator— was shown to be consistent with 
an excess quasiparticle population inferred from noise 
measurements^. There is a vivid debate on the ques- 
tion of the origin of those excess quasiparticlesii^i^, which 
mainly focuses on reducing the influence of the environ- 
ment on the devices under study. Here we show for su- 
perconducting aluminium resonators that the environ- 
ment is well enough under control in our experimental 
setup to reveal a new source of quasiparticles, namely 
the microwave readout power of these devices. We show 
that the saturation in the number of quasiparticles at 
low temperature (100-150 mK), as inferred from noise 
measurements, decreases from 72 to 17 ^Tn~^ with a 6 
dB decrease of the microwave power. The quasiparticle 
lifetime increases accordingly from 1.4 to 3.5 ms. 

Microwave resonators are popular devices in radiation 
detection^ and circuit quantum electrodynamics''. The 
two quadratures of the microwave field in such a res- 
onator are proportional to the real and imaginary part 
of the conductivity of the superconductor ai — ia2- The 
real part corresponds to dissipation in the quasiparticle 
system and the imaginary part to the kinetic inductance 
of the condensate^. We have shown recently that the 
real part of the conductivity shows quasiparticle num- 
ber fluctuations'*. When two quasiparticles recombine, a 
Cooper pair is formed and when a Cooper pair is broken 
it leaves two quasiparticles. Therefore, the superconduct- 
ing condensate fluctuates as well, and one would expect 



to see these fluctuations in the reactive response of the 
microwave resonator. However, this reactive response 
is obscured by the response of two-level fluctuators in 
the dielectrics surrounding the resonator^. Therefore, we 
study here the correlation between the dissipative and 
reactive part of the conductivity in an aluminium res- 
onator. We observe correlated fluctuations in the dissi- 
pative and reactive parts of the response, which proves 
the correlated nature of fluctuations in the quasiparticle 
system and the condensate. The correlation results in a 
measurement of fluctuations in the superconductor down 
to the vacuum noise level, even with a conventional am- 
plifier. The number of quasiparticles and their lifetime 
are extracted from the fluctuation spectra. 

In thermal equilibrium, the average number of quasi- 
particles per unit volume in a superconductor fol- 
lows an exponential temperature dependence: Uqp cx 
y/AltBT exp{—A/kBT). The average quasiparticle life- 
time has the inverse temperature dependence^: Tqp oc 
y^l/ AksT exp{ A/ ksT). A is the energy gap of the su- 
perconductor, T the temperature and ks Boltzmann's 
constant. Two quasiparticles with opposite spins and 
momenta can be generated from a Cooper pair by a 
phonon with an energy larger than 2 A. When two quasi- 
particles recombine into a Cooper pair, a phonon is emit- 
ted. These processes are random processes in equilibrium 
and lead to fluctuations of the number of quasiparticles 
around the average. The power spectral density of these 
fluctuations shows a Lorentzian spectrum, given hy^ 
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with / the frequency and Nq. 



oV, with V the vol- 



Electronic mail: |p.j.devisser@tudelft.nl| 



ume of the system. Since the temperature dependences 
of Nqp and Tqp are exactly opposite, NqpTqp is constant 
over temperature. This formulation was verified through 
fluctuations in the quasiparticle current of a Cooper pair 
bo»ii. Recently we have observed these fluctuations in 
the dissipative response of a microwave resonator—. As 
discussed above, equilibrium fluctuations in the number 
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of quasiparticles should coincide with fluctuations in the 
condensate, which should show up in the reactive re- 
sponse of the resonator. 

To measure the complex conductivity, a 40 nm thick Al 
film was sputter-deposited onto a c-plane sapphire sub- 
strate. The critical temperature is 1.11 K, from which 
the energy gap A = l.TBfcsTc = 168 fieV. The film was 
patterned by wet etching into half wavelength, coplanar 
waveguide resonators. The resonator has a central strip 
volume of 1.0 • 10^ /zm'^ and resonates at 6.62 GHz. The 
sample is cooled in a pulse tube pre-cooled adiabatic de- 
magnetization refrigerator to 100 mK, with a box-in-box 
configuration and coax cable filters for thorough stray 
light shielding, crucial for these measurements. More de- 
tails on the setup are given in Ref. [H- 

The complex transmission of the microwave circuit is 
measured with a quadrature mixer as a function of fre- 
quency and traces out a circle in the complex plane. The 
resonator amplitude. A, measured from the circle center, 
is proportional to ai and therefore called the dissipation 
quadrature. The phase, 9, is proportional to (72 and is 
also called the frequency quadrature. The responsivi- 
ties of amplitude and phase to a change in the number 
of quasiparticles are determined experimentally as de- 
scribed in Ref. [H, which leads to dA/dNqp = — 5 x 10^^ 
and de/dNqp = 4 x 10"^ at 100 mK. At the end of this 
Letter we will discuss the reliability of this method for 
a readout-power dependent quasiparticle density. The 
cross power spectral density due to correlated quasipar- 
ticle number fiuctuations in the resonator amplitude and 
phase is given by 



dA-d9/dN^p 

1 + {2TTfTresY 



(2) 



which is only different from the amplitude or phase power 
spectral density by the responsivity factor, which would 
be {dA/ dNqpf and [dO/dNqpf for the amplitude and 
phase spectra respectively. Tres is the resonator ringtime 
given by Tres = ~ 2 /zs. Because the amplitude re- 
sponsivity to quasiparticles is negative, we expect that 
the correlation of the quasiparticle fiuctuations in ampli- 
tude and phase is negative. 

We have measured the fluctuations in the resonator 
amplitude and phase as a function of time at the res- 
onant frequency. Occasionally peaks occur in the time 
domain data due to high energy impacts, which are fll- 
tered out of the spectral analysis as discussed in Ref. 0- 
The power spectral densities of amplitude and phase are 
calculated by taking the Fourier transform of the auto- 
correlation of the time domain signals. The cross power 
spectral density is calculated by Fourier transforming the 
cross-correlation function of amplitude and phase. If the 
direction of the fluctuations with respect to the resonant 
circle in the complex plane is offset, one may convert 
phase noise into the amplitude direction. By simulating 
different orientations, we estimate the statistical error 
in the orientation to be ±0.28°, which leads to an uncer- 
tainty of ±0.7 dB in the level of the cross power spectum. 
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FIG. 1. (Color online) Amplitude, phase and cross power 
spectral densities of the resonator as a function of frequency 
at 120 mK and a microwave power of —75 dBm. Because of 
the log-scale the positive and negative parts of the cross power 
spectral density are plotted separately. The dashed line is a 
single-timescale Lorentzian fit. 



The amplitude, phase and cross power spectral densi- 
ties at 120 mK are shown as a function of frequency in 
Fig. [TJ We first compare the different spectra. To start 
with, the phase noise is orders of magnitude higher than 
the amplitude noise, which is due to the response of two- 
level-system (TLS) fiuctuators in the dielectrics to the 
electric field?-. Therefore observing quasiparticle fluctu- 
ations in the phase-only spectrum is nearly impossible. 
The fiat level of — 106 dBc/Hz is due to the amplifier noise 
and attenuation in between the sample and the (HEMT) 
amplifier and corresponds to a system noise temperature 
of 7 K. It was recently shown for similar microwave res- 
onators that there is no TLS noise in the amplitude down 
to the vacuum noise^*. Therefore, if the system noise is 
subtracted, the quasiparticle signature becomes visible 
in the amplitude spectrum^. The cross power spectrum 
shows no TLS noise or amplifier noise, which shows that 
these contributions are uncorrelated. 

We now look closer at the cross power spectral den- 
sity. Part of the cross power spectrum (up to 10 kHz) 
is negative (blue) as expected for quasiparticle fiuctua- 
tions. This part of the spectrum is real, meaning that 
the quasiparticle fluctuations enter amplitude and phase 
without relative delay. A small part at higher frequency 
has equal real and imaginary parts, of which the real part 
is positive (red). The negative part consists of two roU- 
offs. The first roll-off is at the quasiparticle lifetime (r 
2 ms, / « 80 Hz) as shown by the dashed line in Fig. 
[T] As a function of temperature, the lifetime from the 
cross spectra is the same as from amplitude-only spectra 
(Ref. W), consistent with the framework of quasiparticle 
number fluctuations. The difference in the level of the 
cross and amplitude spectra (10 ± 1) is due to the dif- 
ference in amplitude and phase responsivity. The first 
conclusion of this Letter is that the correlated noise in 
the amplitude and phase of the resonator is due to cor- 
related fiuctuations in the quasiparticle system and the 
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FIG. 2. (Color online) Cross power spectral density of the 
resonator amplitude and phase as a function of frequency at 
two different microwave powers and temperatures of 120 and 
200 mK. The dashed lines are single-timescale Lorentzian fits. 
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FIG. 3. (Color online) (a) Number of quasiparticles and the 
quasiparticle lifetime as a function of temperature for four 
different microwave readout powers, (b) Quasiparticle den- 
sity and quasiparticle lifetime as a function of readout power. 
Each point is a weighted average of the values from 100-150 
mK as shown in (a). 



superconducting condensate. 

The second roll-off in Fig. [T] is at a shorter timescale 
(r w 100 /iS, / sa 1.5 kHz) and has a much lower (a 
factor 25) noise level. We interpret this second roll-off 
in the spectrum as a signature of phonon fluctuations. 
This phenomenon requires a more extensive discussion 
which we will publish separately. The positive part of 
the spectrum is small and only visible close to the res- 
onator ring-time (which also determines the roll-off fre- 
quency of the phase spectrum), which we attribute to 
phase-amplitude mixing due to slight detuning from the 
resonant frequency during the measurement^^. The sign 
of this contribution varies between different resonators, 
where the quasiparticle contribution is always negative. 

We observe that around 10 kHz the spectral density 
drops to below —120 dBc/Hz, which corresponds to the 
vacuum noise \hf /kB- Note that the vacuum noise is 
not a physical hmit here, but only used for comparison. 
Thus correlating amplitude and phase means a factor of 
25 improvement with respect to the amplifier noise level 
when measuring amplitude or phase only, proving the 
high sensitivity of this method in measuring quasiparti- 
cle fluctuations. For microwave resonators used as pho- 
ton detectors and limited by uncorrelated noise, we en- 
vision an improvement in sensitivity if one reads out the 
detector by correlating the amplitude and phase signal. 

In Figure [2] cross spectra are shown at 120 mK and 
200 mK for two different microwave readout powers. We 
observe that at 120 mK the roll-off frequency increases 
with increasing microwave power. This behaviour is ob- 
served up to 190 mK. At 200 mK, the spectra at different 
powers are similar, which marks the point where thermal 
quasiparticles start to dominate. We extract the number 
of quasiparticles from the cross power spectral density by 
using Eqs. [T] and [2] and the quasiparticle lifetimes as de- 
termined from the roll-off frequency of the spectra. The 
number of quasiparticles and the quasiparticle lifetime 
are plotted for four different readout powers as a func- 



tion of temperature in Fig. The error bars represent 
statistical errors as obtained from the fits and, for iV^p, 
the orientation uncertainty described above. The level at 
which the number of quasiparticles saturates clearly de- 
creases with decreasing readout power. The quasiparticle 
lifetime increases with decreasing power, consistent with 
the decreasing number of quasiparticles. To get a better 
estimate of the saturation levels, we average Uqp and Tqp 
from 100-150 mK and plot the averages as a function of 
power in Fig. ^p. We conclude that the microwave read- 
out signal is the main source of excess quasiparticles at 
low temperature for superconducting resonators. 

In the broader field of superconducting quantum 
circuitsi^ excess quasiparticles are mainly attributed 
to environmental effects, which we strongly reduce by 
shielding our sample box and filtering the cable&iS,. The 
fact that we reveal a microwave power dependence of the 
number of quasiparticles proves that our setup is light- 
tight to at least the lowest measured quasiparticle den- 
sity (less than 0.1 fW of stray-light^). We note that 
since qubits are operated in the limit of a few microwave 
photons, it is unlikely that the excess quasiparticles in 
these systems are also due to the microwave power. 

In the simplest picture, the power that is needed to 
create a certain number of quasiparticles is given by 
Pqp — NqpA/Tqp. Bascd on Pqp and an estimate of 
the power that is dissipated in the quasiparticle sys- 
tem, Pdiss, we can ascribe an efficiency r]read to the pro- 
cess of quasiparticle creation due to the microwave read- 
out power given by Ty^eod = Pqp/Pdiss- Pdiss is given 
by Pdiss = XcXqpPreadl^, with Pread the applied mi- 
crowave power {not the internal power in the resonator). 

Xc = -^Q. is the coupling efficiency and Xqv = Qi/Qi,qp 
is the fraction of dissipated power that goes into the 
quasiparticle system. For this resonator the coupling 
quality factor is Qc = 3.9 x 10^ and the internal qual- 



4 



ity factor is Qi = 1.65 x 10^ at 100 mK and, which 
hardly change with power. Since we can measure the 
number of quasiparticles, we can estimate the quahty 
factor due to quasiparticles, Qi^qp, for each power, which 
ranges from 6 x lO'^ at Pread = —73 dBm to 4 x 10^ at 
Pread — ~79 dBm. Givcn these numbers we calculate 
Vread = 3.3 ± 1.3 X 10~*. The sensitivity of microwave 
resonators used as kinetic inductance detectors is usually 
expressed by the noise equivalent power (NEP). The 
NEP due to quasiparticle number fluctuations^*' can be 
expressed as NEP = "^y/Nqp/r^p = ^^/VreadPdzss^, 
with « 0.6 a conversion efficiency of optical energy 
into quasiparticles. The measured value of rjread is pleas- 
ingly low in this context. At the lowest measured readout 
power we get NEP = 2x 10"^^ W/Hz-i/^. 

So far we have not touched upon the mechanism with 
which the readout signal leads to excess quasiparticles. 
We may explain the excess quasiparticles by Joule heat- 
ing due to the microwave power, leading to an elevated 
steady state temperature of the quasiparticle system. To 
that end, we use the model described in Ref. IT3, in which 
heat transport is assumed to be limited by electron- 
phonon coupling. The model parameters are the same as 
the measured parameters of our device^^. We find that 
Vread ranges from 1 x 10"** at the lowest readout power 
of -79 dBm to 9 x lO"* at -73 dBm. Experimentally 
Tiread IS coustaut withiu the uncertainty. Thus the order 
of magnitude of r]read in the simulations agrees with the 
measurements, but the power dependence is different. 

The question remains how quasiparticle creation by 
the microwave field can be understood microscopi- 
cally. A microscopic picture could be a change of 
the quasiparticle distribution function due to microwave 
absorptioni^Ti^i and consequently a change in the com- 
plex conductivity^^'^^ or an altered density of states due 
to the microwave ficld^'^. The redistribution of quasi- 
particles could lead to Cooper pair breaking. The next 
step to unravel the physical mechanism of the microwave 
power dependent quasiparticle density will therefore be a 
comparison of these models to resonator measurements. 
In this context we note that in Fig. [3}3, the product 
NqpTqp is not completely constant as a function of read- 
out power. Our analysis here is based on thermal quasi- 
particles, from which we derive the responsivity to quasi- 
particles. As the quasiparticle distribution may be non- 
thermal, the complex conductivity and therefore the re- 
sponsivity to the number of quasiparticles may change. 
This could affect the derivation of Nqp from the noise 
level. The measurement of Tqp is fairly rigid, since it is 
obtained from the roll-off frequency only. 

In summary, we have measured correlated fluctuations 
in the quasiparticle system and condensate of a super- 
conductor, which show up in the amplitude and phase 
of an aluminium microwave resonator respectively. From 
the correlated noise spectra, we determine the number 
of quasiparticles and their lifetime, which both saturate 
at temperatures below 160 mK. The level of this satu- 
ration is microwave power dependent, showing that the 



microwave readout power leads to excess quasiparticles. 
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